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1 Introduction 



1.1 FRW cosmology 

As indicated by observations, the universe is assumed to be homogeneous and 
isotropic on large scales. The line element is: 



ds 2 = dt 2 -a 2 (t) 



+r 2 (d6 2 + sin 2 6dtf) 

1 — kr z 



(1) 



where the function a(t) is known as the scale factor, and the parameter k — 
— 1,0, +1 labels an open, flat or closed universe. The Einstein equations give: 



d\ k 8nG . . 



where p(t) is the energy density of the cosmological fluid. Conservation of 
energy-momentum yields 

p + 3H(p + p) = (3) 

where H = a/a is the Hubble parameter and p denotes the pressure of the 
cosmological fluid. The properties of the cosmological fluid are given by its 
equation of state: 

P = P(p) (4) 

which must be specified to determine the expansion rate of the universe. If the 
fluid is relativistic (eg. photons): p = p/3; while if it is dust: p — 0; and for 
vacuum energy (cosmological constant): p = —p. 

If the fluid is in thermal equilibrium, the evolution of the temperature of 
the fluid can be derived by using the conservation of entropy which gives 

a(t)T — constant . (5) 

The temperature of the cosmic microwave background radiation (i.e. cosmo- 
logical photons as opposed to those produced by stars and other astrophysical 
sources) at the present epoch is measured to be 2.7°K. 

It is sometimes convenient to work with the "conformal time" r instead of 
the "cosmic time" t. The relation defining r in terms of t is: 

dt = a(t)dr . (6) 

In terms of the conformal time, the FRW line clement is: 



ds 2 = a 2 (t( T )) 
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dr 2 - + r 2 (d9 2 + sin 2 Bd<$ 2 ) 

1 - kr z 



(7) 
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1.2 Cosmological phase transitions 

We consider a field theory of scalar, spinor and vector fields: 

L = L B +L F (8) 

with the bosonic Lagrangian: 

L B = ^D^D^i - V($) - -F%F^ a (9) 

where $i are the components of the scalar fields. The fermionic Lagrangian 
for a fermionic multiplet ^ is: 

Lp = iV^D^ - *r,*$, . (10) 

In addition we have the following definitions: 

D^ = d^-ieA^T a (11) 

where the T a are group generators; 

F; v = d^A a y - d v A« + +ef abc AlAl (12) 

where, A® are the gauge fields. 

If the expectation values of the scalar fields are denoted by $oi, then the 
mass matrices of the various fields are written as: 



2 d 2 V 
My = 



, scalar fields (13) 

<3>— "l»o 



m = Ti&oi , spinor fields (14) 
where, the I\ are the Yukawa coupling matrices, and 

M lb = e 2 (T a T b ) t ^ ^ 03 ; , vector fields (15) 

Then the finite temperature, one-loop effective potential is: 

V eff (* ,T) = y($ ) + ^T 2 - ^AAT 4 (16) 

where 

N = N B + IaTf (17) 



is the number of bosonic and fermionic spin states, and 

M 2 = Tr/i 2 + 3TrM 2 + ^Tr (^m^m) . (18) 

Note that M 2 depends on the expectation value <&o through the defining equa- 
tions for the mass matrices given above. Then, for example, M 2 will contain 
a term proportional to Tt($q). 

For us the important feature of the effective potential is that it can lead 
to cosmological phase transitions. If there are scalar fields with negative mass 
squared terms in V(&), the contributions from the Ai 2 T 2 term in the effective 
potential can make the effective mass squared positive for these fields if the 
temperature is high enough. Therefore when the universe is at a high temper- 
ature, the effective squared mass is positive and the minimum of the potential 
is at $ = 0. In the more recent cooler universe the effective mass squared is 
negative and the minima of the effective potential will occur at non-zero values 
of <!>o- That is, the scalar fields will acquire vacuum expectation values. This 
is the phenomenon of spontaneous symmetry breaking and, if it occurs, will 
manifest itself as a cosmological phase transition. 

1.3 Vacuum manifold 

In general the field theoretic action under consideration will be invariant under 
some transformations of the fields. The set of all such symmetry transforma- 
tions form a group which we denote by G. 

Next consider the situation where a scalar field acquires a vacuum expec- 
tation value (VEV) <j>o- Then all transformations of (j>o by elements of G will 
also be legitimate VEVs of the scalar field. However, not all elements of G will 
have a non-trivial effect on <j) . There will be a subgroup of G which will have 
a trivial action on (j) . This is the unbroken subgroup of G and we denote it 
by if. 

Now we wish to determine the set of possible VEVs of the scalar field. 
This set is given by the elements of G that have a non-trivial action on (j)Q. 
Now consider an element of g £ G that has a non-trivial action on </> and 
transforms it to (pi. Then every element of the form gh where h £ H, also 
takes 4>o to 4>i ■ Therefore the non-trivial transformations are given by the left 
cosets of H, and the space of all non-trivial transformations is called a coset 
space and denoted by G/H. Therefore the vacuum manifold of the theory is 
G/H. 

As we shall see, topological defects occur due to non-trivial topology of 
the vacuum manifold. 
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Note that the topology does not care if the symmetries are local or global. 
For most of these lectures, I will consider local (gauged) symmetries. 

2 Domain walls 

2.1 Topology: ttq 

Domain walls occur when the vacuum manifold has two or more disconnected 
components. For example, G can be a discrete group like Z 2 which can be 
broken completely (H = 1). Then G/H consists of just two points. 

The topology of various manifolds has been studied by mathematicians 
and they characterize the type of topology by stating the homotopy groups of 
the manifold. The n th homotopy group of the manifold G/H is denoted by 
ir n (G/H) and denotes the group that is formed by considering the mapping of 
n— dimensional spheres into the manifold G/H. Each clement of the group ir n 
is the set of mappings that can be continuously deformed into one another. 

In the case when the vacuum manifold has disconnected components, 
tt (G/H) is non-trivial since there are points (zero dimensional spheres) that 
lie in different components that cannot be continuously deformed into one 
another. Therefore domain walls occur whenever ttq(G/H) is non-trivial. 

2.2 Example: Z 2 

Consider the Z 2 Lagrangian in 1+1 dimensions labeled by (t, z) 

L= 1 -(d^r-^- v r (19) 

where <j>(t, z) is a real scalar field - also called the order parameter. The La- 
grangian is invariant under 0^—0 and hence possesses a Z2 symmetry. For 
this reason, the potential has two minima: <f> = ^V, an d the "vacuum manifold" 
has two-fold degeneracy. 

Consider the possibility that (p = +77 at z = +00 and <f> = —77 at z = —00. 
In this case, the continuous function <f>{z) has to go from — 77 to +77 as z is 
taken from —00 to +00 and so must necessarily pass through = 0. But 
then there is energy in this field configuration since the potential is non-zero 
when = 0. Also, this configuration cannot relax to either of the two vacuum 
configurations, say (j>(z) = +7], since that involves changing the field over an 
infinite volume from — 77 to +77, which would cost an infinite amount of energy. 

Another way to see this is to notice the presence of a conserved current: 

f = e^d v 4> 



6 



where fj,, v = 0, 1 and e^" is the antisymmetric symbol in 2 dimensions. Clearly 
j M is conserved and so we have a conserved charge in the model: 



Q = I dzf = 0(+oo) - <f>(-ac) 



For the vacuum Q = and for the configuration described above Q = 1. So 
the configuration cannot relax into the vacuum - it is in a different topological 
sector. 

To get the field configuration with the boundary conditions 0(±oo) = ±rj, 
one would have to solve the field equation resulting from the Lagrangian ([l9|). 
This would be a second order differential .equation. Instead, one can use the 
clever method first derived by BogomolnyiEJ and obtain a first order differential 
equation. The method uses the energy functional: 

E = J ' dz[\{d t ct>) 2 + l(d z <f>) 2 +V(4>)] 

dz[\{d t <p) 2 + \{d z <t> - Vm^) ) 2 + JW($jd z ci>] 



4>{+oc) 



dz[-(d t 4>) 2 + -(d z 4> - y/2V&) ) 2 ] + / dep'^WW) 



Then, for fixed values of 4> at ±oo, the energy is minimized if 

d t (f> = 

and 

d z <t> - ^/2V{cf>) = . 

Furthermore, the minimum value of the energy is: 

( + oo) 



E mm = I d^'^fWW) 

In our case, 



which can be inserted in the above equations to get the "kink" solution: 



= rytanh^y^r?^ 



for which the energy per unit area is: 



2v2 i— , 2\/2 



•— \A//' —^-^ (20) 



where m = VA?7 is the mass scale in the model (see eq. ([l9|)). Note that the 
energy density is localized in the region where <p is not in the vacuum, i.e. in 
a region of thickness ~ m" 1 around z = 0. 

We can extend the model in eq. ( |l9| ) to 3+1 dimensions and consider the 
case when <f> only depends on z but not on x and y. We can still obtain the 
kink solution for every value of x and y and so the kink solution will describe 
a "domain wall" in the xy— plane. 

At the center of the kink, </> = 0, and hence the Z2 symmetry is restored in 
the core of the kink. In this sense, the kink is a "relic" of the symmetric phase 
of the system. If kinks were present in the universe today, their interiors would 
give us a glimpse of what the universe was like prior to the phase transition. 



2.3 SU(5) 

An example that is more relevant to cosmology is motivated by Grand Unifi- 
cation. Here we will consider the SU(5) model: 

L = Tr(^$) 2 - ^Tr(X^X^) - V($) (21) 

where, in terms of components, $ = <I> a T a is an SU(5) adjoint, the gauge field 
strengths are X^ v = X^ v T a and the SU(5) generators T a are normalized such 
that Tv(T a T b ) = 6 ab /2. The definition of the covariant derivative is: 

= - ie[X„, $] (22) 

and the potential is the most general quartic in $: 

V($) = -m 2 Tr($ 2 ) + /i[Tr($ 2 )] 2 + ATr($ 4 ) + 7 Tr($ 3 ) - V , (23) 



where, Vo is a constant that we will choose so as to set the minimum value of 
the potential to zero. 

The SU(5) symmetry is broken to [SU(3) x SU(2) x U(l)]/Z 6 if the Higgs 
acquires a VEV equal to 

$o = ^=diag(2,2,2,-3,-3) (24) 
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where 



For the potential to have its global minimum at $ = $o, the parameters are 
constrained to satisfy: 

A > , A' > . (26) 
For the global minimum to have V($o) = 0, m eq. (|2^ ) we set 

Vo = -j V 4 - (27) 

The model in eq. ( ^l|) does not have any topological domain walls because 
there are no broken discrete symmetries. In particular, the Zi symmetry under 
$ — > — $ is absent due to the cubic term. However if 7 is small, there are walls 
connecting the two vacua related by $ — * — $ that are almost topological. 
In our analysis we will set 7 = 0, in which case the symmetry of the model 
is SU(5) x Z2 and an expectation of <E> breaks the Z2 symmetry leading to 
topological domain walls. 

The kink solution is the Z2 kink along the <E>o direction (see eq. (pi])). 
Therefore: 

$ fc = tanh(crz)$ (28) 

with a = m/y/2 (see eq. (pSj)), and all the gauge fields vanish. It is straight- 
forward to check that $fc solves the equations of motion with the boundary 
conditions $(z = ±00) = ±$o- 

The energy per unit area of the kink is (see eq. (p0|)): 

2V2m 3 , , 

- — -y ■ (29) 

The existence of a static solution to the equations of motion only guar- 
antees that it is an extremum of the energy but this extremum may not be 
a minimum. To determine if the kink is a minimum energy solution we need 
to examine its stability under arbitrary perturbations. (As far as I know, a 
Bogomolnyi type analysis has not been constructed for the SU(5) model.) 

Here we will examine the stability of the kink under general perturbations. 
So we write: 

$ = + * (30) 

Since the kink solution is invariant under translations and rotations in the 
xy— plane, it is easy to show that the perturbations that might cause an in- 
stability arise from perturbations of the scalar field and can only depend on z. 
Therefore we may set the gauge fields to zero and take * = ^(t, z). 
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The Z2 kink is stable and hence we can restrict the scalar perturbations to 
be orthogonal to. . The perturbative stability .of the kink has been studied 
by Dvali et. al. □, and, Pogosian and Vachaspatiu. Here we only consider the 
off-diagonal perturbation 

* = i)T = ^diagfV 1 , 0,0,0) , (31) 

where r 1 is the first Pauli spin matrix. 

Next we analyze the linearized Schrodinger equation for small excitations 
ifi = 'tJjo(z)exp(~iLut) in the background of the kink: 

- m 2 + 4>l{z){h + Ar)]^ = ^Vo , (32) 

where (f>k = tanh(crz) and r = 7/30. The kink is unstable if there is a solution 
to eq. (E32I) with a negative u 2 . Substituting eq. J2§|) into eq. ((32|) yields: 



[-d 2 + m 2 (tanh 2 (az)-l)}i[>o=Lj 2 Tp . (33) 

This equation has a bound state solution ipQ oc sech(crz) with the eigenvalue 
lo 2 = —m 2 /2. Since this result is independent of the parameters in the poten- 
tial, we conclude that the kink in SU(5) is always unstable. 

So we still need to find the topological domain wall solution in the model. 

The domain wall solution is obtained if we choose the gauge fields to vanish 
at infinity and the scalar field to satisfy the boundary conditions: 

Hz = -co) = $- = -^=diag(3, -2, -2, 3, -2) 
2v 15 



and 



$(z = +co) = $+ = — L=diag(2, -3, 2, 2, -3) 

= ^yj(A3 + r 3 ) + |(r-V5A 8 ) . (35) 



Here A3, A§, T3 and Y are the diagonal generators of SU(5): 



A 3 = idiag(l,-l,0,0,0) , (36) 
' diag(l, 1,-2, 0,0) , (37) 



2^3 
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r 3 = -diag(0, 0,0, 1,-1) , 



(38) 



F=— J=diag(2,2,2,-3,-3) . (39) 
2y 15 

Alternately, we could have chosen the boundary conditions to be like those of 
the kink with $(z = +oo) = — $(z = — oo). However, then the solution for the 
domain wall will not be diagonal at all z. We prefer to use the above boundary 
conditions so that the solution is diagonal throughout. 
The domain wall solution can be written as 

$ DW (z) = a(z)X 3 + b(z)X 8 + c(z)t 3 + d{z)Y . (40) 

The functions a, b, c, and d must satisfy the static equations of motion: 

a" = l-m 2 + (h+— )d 2 + 
5 

+ 



(h + -)(a 2 + b 2 ) + hc 2 ]a 

2Xabd ... 

(41) 



b" 



-m 2 + (h+^)d 2 
5 



(h+-){a 2 + b 2 ) + hc 2 ]b 

7E [a ~ b) 



(42) 



9A 



A, 



+ (h+-)d 2 + (h + -) 
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2' 
h(a 2 



b 2 )]c 



(43) 



d" = [-m 2 + (h+ 7 ^)d 2 + (h+^)(a 2 + b 2 ) 
30 5 

9A, ,,, Ai, 2 6 2 , , AA . 

+ ( h+ w ) ] + ^E {a 3" } ' ( } 

where primes refer to derivatives with respect to z. For reference, the kink so- 
lution (eq. (p8|)) corresponds to a(z) = = b(z) — c(z) and d(z) = i]ta,nh(az). 
The equations of motion for b and c and can be solved quite easily: 

b(z) = -Vbd{z) , c(z) = a{z) . (45) 
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This is consistent with the boundary conditions in eqs. (34) and (35). In 
addition, we require 



/ 5 r\ 
alz = ±oo) = +W\ — , d(z = ±oo ) = ± — 
V 12 6 



(46) 



Then the remaining equations can be written in a cleaner form by rescaling: 

(47) 



. / x /12 a „ , , d „ 

A(z) = \ , D(z) = 6~, Z = mz. 

V 5 T) rj 



This leads to 



A" = 



D 



-1 



(1-P) g2 , (4+P) 



-l+p£> 2 + (l-p)A 2 



4 2 



D 



.4 



where primes on A and _D denote differentiation with respect to Z , and 



5A 
12A 7 



(48) 
(49) 

(50) 



Note that p € [1/6, oo) because of the constraints in eq. (|26|). The boundary 
conditions now are: 



A(z = ±oo) = +1 , D(z = ±oo) = ±l 



(51) 



This system of equations has been solved by numerical relaxation and a 
sample solution is shown in Fig. [I]. To find an approximate analytical solution, 
assume that |A"/j4| << 1 is small everywhere. This assumption will be true 
for a certain range of the parameter p which we can later determine. Then the 
square bracket on the right-hand side of eq. (B8h is very small. This gives: 



.4 



4 + p 



1 - ^ —D 2 



-l 1/2 



(52) 



We insert this expression for A in eq. (|49|) and obtain the kink-type differential 
equation: 

D" = q[—l + D 2 ]D , (53) 



where 



6A 



6p — 1 
p + A ~ A + 60A' 



(54) 
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Figure 1: Numerical solution for the domain wall in the case X = 1, h = —0.2 (p = 2.25). 
The solid line shows a(z) and the dashed line shows d(z). 



and the solution is: 



The parameter q lies in the interval [0, 6]. For q = 1 (i.e. p = 1) it is easy to 
check that this analytical solution is exact. 

We can now check that our assumption |^4"/^4| << 1 is self-consistent 
provided p is not much larger than a few. 

The energy density for the fields A and D can be found from the La- 
grangian in eq. (|2|) together with the ansatz in eq. (fit)]), the solution for b 
and c in eq. ( pl5| ) and the rescalings in eq. (|47]). The resulting expression for 
the energy per unit area of the domain wall is: 




(55) 



M DW 



12A' 



m 




dZ[5A' 2 +D' 2 +V(A,D)} 



(56) 



where 



V(A,D) 



5A 2 



+ ^D i + {l-p)A 2 D 2 +3 . 




(57) 
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Table 1: Size distribution of + clusters found by simulations on a cubic lattice. 



Cluster size 


1 


2 


3 


4 


6 


10 


31082 


Number 


462 


84 


14 


13 


1 


1 


1 



The energy can be found numerically. However, here we will find an approx- 
imate analytic result. We can insert the approximate solution given above in 
eq. (^6|) but this leads to an expression that is not transparent. Instead it is 
more useful to consider another approximation for A and D: 

D ~ tanh^y^zj . (58) 

(This approximation is exact for p = 1.) A straightforward evaluation then 
gives: 

Md\V approx = Mfc-— (59) 

where, Mk is given in eq. (p9|). 

It can be shown for a range of parameters that this domain wall solution 
is perturbatively stable. Numerically we find that it is lighter than the kink 
for all values of p. (Eq. ( |59| ) shows it to be lighter for p < 36.) 

The interior of the domain wall has less symmetry than the exterior. So the 
domain wall does not contain a trapped region of the early universe. Instead 
it is a region that may be to our future where the SU(3) symmetry is broken 
down to SU(2) x U(l). In this sense, the domain wall is a relic from the early 
universe but not a relic of the unbroken symmetry of the theory 



2-4 Formation 

The properties of the network of domain walls at formation has been deter- 
mined by numerical simulations. The idea behind the simulations is that the 
vacuum in any correlated region of space is determined at random. Then, if 
there are only two degenerate vacuua (call them + and -), there will be spatial 
regions that will be in the + phase and others in the - phase. The boundaries 
between these regions of different phases is the location of the domain walls. 

By performing numerical simulations, the statistics shown in Table I was 
obtained!! 
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The data shows that there is essentially one giant connected + cluster. By 
symmetry there will be one connected - cluster. In the infinite volume limit, 
these clusters will also be infinite and their surface areas will also be infinite. 
Therefore the topological domain wall formed at the phase transition will be 
infinite. 



2.5 Evolution 

Once the inifinite domain wall forms, it tries to straighten out since it has 
tension - in this way, the wall can minimize its energy. In the cosmological 
setting, the wall will be stretched out by the Hubble expansion. In addition, 
the wall will interact with the ambient matter and suffer friction. Sometimes 
different parts of the walls will collide and interact. The interaction leads to 
the reconnection of colliding walls. 

The details of the evolution of the walls are quite complicated. However, 
the simplest cosmological scenario is easy to analyze and leads to an interesting 
bound. The basic idea is to use causality to get an upper bound on the energy 
density in domain walls at any epoch. Causality tells us that the different 
domains that are separated by walls cannot smooth out faster than the speed 
of light. Therefore every causal horizon must contain at least one domain 
wall at any epoch and so the energy density in domain walls pnw obeys the 
following bound: 

^ energy in one wall in horizon at 2 a 
horizon volume i 3 t 

where a is the energy per unit area of the wall. If we require that the energy 
density in walls be less than the present critical density of the universe p cr = 
3H 2 /8nG with H ~ 70 km/s/Mpc we get the constraint 

a< ^' (61) 

where to ~ 10 17 sees is the present epoch. 

If we use eq. (^0|) to connect a to the symmetry breaking scale and take the 
coupling constant to be order one, we find that the symmetry breaking scale 
is constrained to be less than about a GeV. A stronger bound is obtained by 
realizing that the energy density fluctuations in domain walls should not cause 
fluctuations in the cosmic microwave background radiation (CMBR) greater 
than 1 part in 10 5 . This leads to the constraint that r\ should be less than 
a few MeV. Particle physics at such low energy scales - such as the standard 
model - do not show any phase transitions involving domain walls and so we 

15 



do not expect to have domain walls in the universe. Another way to view the 
constraint is that if a particle physics model predicts domain walls above the 
few MeV scale, cosmology rules out the model. 

2.6 Further complexities 

One can evade the constraints derived in the previous section by introducing 
some new cosmological and particle physics elements. 

For example, the constraint that domain walls have to be formed below 
the few MeV scale does not hold if cosmological inflation follows the phase 
transition in which heavy domain walls were produced. In such a scenario, the 
domain walls are inflated away and our present horizon need not contain any 
domain wall. 

In particle physics model building, one might have a situation where a 
discrete symmetry is broken at some high temperature and then restored at 
some lower temperature. In this case, heavy domain walls would be formed 
at the first phase transition and then would "dissolve" at the second phase 
transition. Then there would be a period in the early universe where domain 
walls would be present but they would not be around today. 

Another situation of some interest is when the discrete vacuua are not 
exactly degenerate. In the SU(5) example discussed above, if 7 (eq. (53)) is 
very small, domain walls will be formed at the phase transition. Suppose the + 
phase has slightly less energy density than the - phase. Then the domain walls 
will experience a pressure on them that will drive them towards the higher 
energy - phase. However, if 7 is small, this pressure difference is also small and 
is negligible for the early evolution of the wall system. Only after the walls 
have straightened out to an extent that the curvature forces are less than the 
pressure, will the pressure start to drive the walls and eventually eliminate 
them. 

3 Strings 

3.1 Topology: tti 

If the vacuum manifold has one dimensional closed paths that cannot be con- 
tracted, there are topological string solutions in the field theory. The homotopy 
group iri(G/H) is the group formed by the equivalence classes of paths that 
can be deformed into each other and the group operation joins two paths to get 
another path. Each element of tti(G/H) labels a topologically distinct string 
solution. 
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An example of a field theory in which string solutions exist is based on a 
U(l) (global or local) group which is broken down completely. The coset space 
U(l)/1 is a one dimensional sphere (circle). Closed paths that wrap around 
the circle cannot be contracted to a point and this implies that the U(l) model 
contains string solutions. 

The U(l) model is very important as it is relevant to the simplest su- 
pcrfiuids and superconductors. However, in more complicated systems the 
symmetry groups are larger and calculating the homotopy group can be more 
involved. Fortunately there is a result that simplifies matters immensely for 
most (though not all!) particle physics applications. This is the result that if 
7r„(G) and 7r„_i(G) are both trivial then, 

n n (G/H)=n n - 1 (H) . (62) 

With n = 1 this gives 

m(G/H)=Tro(H) (63) 

provided tt\{G) = 1 = n (G). So if G docs not contain any incontractable 
closed paths and does not have disconnected pieces, then the topologically 
distinct incontractable paths of G/H are given by the disconnected pieces of 
H. So one can usually tell that there are strings in a particle physics model 
simply by checking if the symmetry breaking involves a broken U(l) or if the 
unbroken group contains a discrete symmetry. This result cannot be applied 
to the electroweak model since there we have G = [SU(2) x U(l)]/Z2 and 
7ri(G) = Z^l. 

3.2 Example: U(l) local 
Consider the Lagrangian 

where $ is a complex field and 

Dp — ^ ieAfi 

where, is an Abelian gauge potential. This model has a U(l) gauge sym- 
metry since it is invariant under 

$ -> e* e $ , — > + -d„6 (64) 

e 

Now since 7Ti(Z7(l)/l) = Z, the model has string solutions labeled by an integer 
(the winding number). 
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The string solutions correspond to the non-trivial windings of a circle at 
spatial infinity on to the vacuum manifold. Therefore the boundary conditions 
that will yield the string solution will have the form: 



= r\e 



inO 



e z - {0} 



(65) 



The gauge fields at infinity must be such that the covariant derivative vanishes 
there: 

A$ = , at r = oo . (66) 

Let us now construct the string solution using Bogomolnyi's methodS. The 
energy for static configurations in two spatial dimensions is 



E= d 2 x 



\\D V *\ 



with 

The trick is to write this as 
1 



F($) = ^(|$| 2 -r ? 2 ) 2 



E 



d 2 x 



-\D x *±iD v <S>\* 



'ivy 



-;>) 



d 2 xB 7 



(67) 



(68) 



(69) 



in the special case that A = e 2 (so called "critical coupling" ) 0. The Bogomolnyi 
equations are: 

D x $ + iD,,® = 



B z - V2V = 

and the energy per unit length of the string is: 



fl= 2 V 



d 2 xB z = -r] 2 <f> dOAg = trr] 2 



(70) 
(71) 

(72) 



For non-Bogomolnyi strings, this expression will be multiplied by a factor 
which depends on A/e 2 . Unless this parameter is very large or small, numerical 
evaluations show that the numerical coefficient is of order unity. 

For strings formed at the Grand Unified phase transition r/ ~ 10 16 GeV, wc 
find /i ~ 10 22 gms/cm which is like the mass of a mountain range packed into 
a centimeter of string. When calculating the gravitational effects of strings 
\x always appears in the dimensionless combination G[x. For GUT strings, 
Gfi — 10~ 6 . 



a Unlike in the case of the Zi domain 
the form of the potential given in eq. A 



ill, the Bogomolnyi trick over here works only for 
ij) and with critical coupling. 
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3.3 Semilocal and electroweak strings 

The standard model symmetry breaking is SU (2)l x U(1)y —> U (l) em and does 
not give topological strings. However, it still contains embedded strings which 
can be perturbatively stable for a range of parameters which are not realized 
in Nature. In these lectures I will only briefly describe semilocal stringsd. 

Consider the generalization of the Abelian Higgs model in which the com- 
plex scalar field is replaced by an ST/ (2) doublet <& T = (<j>i,<j>2). The action 



S = d*x 



(73) 



where is the U(l) gauge potential and Y^ v — d^Y v — d v Y^ its field strength. 

This model has symmetry under G = [SU(2) g i b a i x U(l)i oca i\/ Z 2 . Ele- 
ments of the SU (2) gi b a i act on the Higgs doublet, while elements of U{l)i oca i 
multiply the doublet by an overall phase and transform the gauge field in the 
usual manner. The action of the center of SU(2) g i b a i on the Higgs doublet 
is identical to a U(l) g i bai phase rotation by n. This is the reason for the Z 2 



identification. The model in eq. (73) is just the scalar sector of the Glashow- 
Salam- Weinberg (GSW) electroweak model with sin 2 9 W = 1. 

Once $ acquires a vacuum expectation value, the symmetry breaks down 
to H = U(l) global, a s in the GSW model. The vacuum manifold is S 3 . This 
may also be seen by minimizing the potential explicitly. The vacuum manifold 
is described by = n 2 /2 which is a three sphere. Hence there are no 

incontractable paths on the vacuum manifold. Yet, it is possible to perform a 
Bogomolnyi type analysis to find that the configuration: 

* = f(py S ( 1 ) , Y 9 = V M (74) 



qp 



is a solution in cylindrical coordinates (p,0). The profile functions f(fl) and 
v(fl) can be determined by solving the Bogomolnyi equations. 

It is also possible to show that the solution is perturbatively stable when 
2 A < q 2 , neutrally stable at critical coupling (2A = q 2 ) and unstable for 2 A > 
q 2 . 

For the case when the SU(2) is gauged, the string solution continues to 
exist and can be stable for sin 2 9 W > 0.9 if the scalar mass is small. In the 
standard model, however, sin 2 8 W = 0.23 and the string solution is certainly 
unstable. Certain external conditions (eg. a magnetic field) can stablize the 
string. 
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3.4 Formation and evolution 



The formation of strings has been studied by numerical methods. Here one 
throws down U(l) phases on a lattice and checks for topological winding around 
the plaquettes. If the winding around a plaquette is non-trivial, it means that 
a string is passing through it. Then the strings are connected and in this way 
one obtains a network of closed and open strings. 

The lattice simulations show several interesting properties of the initial 
network of strings. First and foremost is the result that most (over 80%) of 
the string density is in strings that are infinite. Secondly, the string loops 
occur in a scale invariant distribution - that is, the number density of loops 
having size between R and R + dR is proportional to dR/R 4 , just as would 
be expected based on dimensional analysis. Thirdly, the long strings have 
a Brownian shape and so the length I of string in a distance R is given by: 



There has been much effort in extending the lattice simulations to make 
them more realistic and closer to what might actually happen at a phase tran- 
sition. However, the essential results of the lattice simulations have held up 
very well. There has also been a lot of theoretical and experimental effort de- 
voted to determining the density of defects after the phase transition. (In the 
simulations the density of defects is set by the lattice spacing.) The statistical 
properties of the network described above are expected to be insensitive to 
the initial string density. The exact fraction of strings in infinite strings might 
vary somewhat. 

Once the string network has formed, it evolves under the forces that we 
discussed in the domain wall case: (i) tension, (ii) friction with ambient mat- 
ter, (iii) Hubble expansion, (iv) intercommuting (reconnection) when strings 
collide, and (iv) energy loss mechanisms such as gravitational wave emission. 
Of these factors, frictional forces are important for a short period (for GUT 
strings) since the Hubble expansion dilutes the ambient matter. The energy 
loss mechanisms are important but to calculate their effect on the string dy- 
namics is a nasty back-reaction problem. The "zeroth order" factors are the 
tension, Hubble expansion and intercommutings. 

Ignoring friction, energy loss and intercommutings for now, the strings 
move according to the Nambu-Goto action in a background FRW spacetime: 



where is the determinant of the two dimensional string world-sheet metric. 
If the string worldsheet is written as: x^(a a ) with a = 0, 1 where (a , a 1 ) are 



I oc R 2 . 




(75) 
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the worldsheet coordinates, then the worldsheet metric is: 

9 { a l } = 9„ud a x»d b x v (76) 

where, g M „ is the four dimensional spacetime metric in which the strings move. 

The freedom of choosing coordinates on the worldsheet (reparametrization 
invariance) can be used to simplify the equations of motion for the string 
derived from eq. (|75|). We choose 

a°=r, ^ = C (77) 

where £ is a parameter along the string at any given instant of the conformal 
time t (see eq. @). Then the equation of motion is: 

x + 2H(l -k 2 )± = e- 1 ^-^ (78) 

where overdots refer to derivatives with respect to r and primes with respect 
to C- 

The equation of motion for a single string allows numerical evolution of 
the network of strings but does not take intercommutings into account. In 
practice, the network is numerically evolved until two strings intersect and 
then the reconnection is done by hand. 

To study the evolution of cosmic strings by solving the equations of motion 
is similar to studying a box of gas by solving Newton's equations of motion for 
each particle in the box and taking care of collisions by hand. For almost all 
cosmological purposes, only the statistical properties of the string network are 
relevant. So one should define "statistical variables" and find their evolution. 
The one scale model for string evolution is precisely such an attempt. 



3.5 One scale model 

As the name suggests, the "one scale model" is based on the assumption that 
there is rPflly one length scale that characterizes the infinite strings in the 
network Em. (The short loops are considered separately.) If we denote this 
scale by L then the distance between strings and the correlation length of a 
single string (the distance out to which it is straight) are both given by L. 
With time L can change and hence we need to derive an equation for L(t). 

The usual way of defining the length L is by considering the total energy 
E in strings lying inside a volume V . Then 
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defines L. The energy changes with time since the Hubble expansion dilutes 
and stretches the string network. Also, the velocity of a string rcdshifts due 
to the Hubble expansion. Long strings also lose energy by self-intersecting to 
form short loops. This gives 

p = -2H(1 + v 2 rms )p - (80) 



where 



rx 



2 



Vims = = ( 81 ) 

is the average rms velocity of strings and c is a parameter that accounts for the 
frequency of self-intersection (the so called "chopping efficiency"). Equation 
( pp| ) can also be written as an equation for L: 

L = HL(\ + vl ms ) + ^ (82) 

(Since the chopping term must vanish if the velocity of the strings vanishes, 
one frequently writes c = cv rms .) 

The network is said to scale if, at every epoch, the network has the same 
statistical properties. Hence the scale L should be a fixed fraction of the horizon 
size t. In this case, the energy density p in long strings will fall off like 1/t 2 , 
exactly as the total matter density does in a radiation or matter dominated 
universe. In a universe with a cosmological constant, scaling will hold in a 
trivial sense since the rapid expansion will dilute the strings to vanishing energy 
density. Scaling will not hold in the epochs during the period in which the 
universe changes from being radiation dominated to matter dominated, or from 
matter dominated to cosmological constant dominated. Both these transitions 
are relevant for examining the cosmological signatures of strings (eg. the cosmic 
microwave background anisotropics), making the problem quite hard. 

To check for scaling we write 

L = 7 (t)t , (83) 
insert this relation into eq. (p0|) to get 



1 = 1 
7 t 



{l-Ht(l + v 2 rms )} 

2 7 



(84) 



There is a fixed point solution (7 = 0) at the point: 



c 1 

7= 21-(Ht)(l + v 2 ms ) • (85) 
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Note that for a power law expansion Ht does not depend on time. The rms 
velocity could depend on time though we can prove that v^ ms = 1/2 in flat 
spacetime suggesting that it may be a reasonable assumption to take it to be 
roughly constant even in the expanding universe. However, we can derive an 
equation for v rms by differentiating eq. (|8l]) with respect to time and then 
using the equation of motion. This gives 



(1-v, 



2 ) 

/•///.- J 



k 



(86) 



where k is a parameter (the "momentum parameter") related to the typical 
radius of curvature of the strings. The second term gives the redshifting of the 
velocity due to Hubble expansion while the first describes the change in the 
velocity due to the tension in the string. 

Equation (|2|) and fl86| ) are the one scale model equations for the string 
network. These equations depend on the chopping efficiency c (or c) and the 
momentum parameter k that need to be specified. The practice has. been to 
use the values found in numerical simulations. Martins and ShellardQ give 

c = 0.23 ±0.04 (87) 

and a velocity dependent value of k 

. 2V2 1-8v 6 
y ' n l + 8v 6 y ' 



4 Cosmological constraints and signatures 

Numerical evidence suggests that the density of strings scales with the expan- 
sion of the universe and so we have: p ~ /i/£ 2 - Compared to the critical density 
in a radiation- or matter-dominated universe p c ~ 1/Gt 2 , the string density 
is at least a factor ~ G/i smaller at all times. So the strings never come to 
dominate the universe. 



4.1 CMB and P(k) 

The energy-momentum of the network of strings causes perturbations in the 
metric which can introduce fluctuations in the cosmic microwave background. 
The metric produced by a straight string is conical with deficit angle 87rG/i and 
photons arriving to us from either side of a moving string have a temperature 
difference of ~ SttG^lv where v is the velocity component of the string trans- 
verse to the direction of propagation of the photons. Therefore, on dimensional 
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grounds, the CMBR anisotropy produced by strings is: 

ST 



(89) 



which for GUT strings is about what is observed. 

Several sophisticated calculations of the CMBR anisotropy induced by 
strings have been done. I would like to draw a distinction between some of 
the analyses that have been done for global strings and others that have been 
done for local strings. Here I will only discuss the specific analysis done by 
Pogosian and me of CMBR anisotropies produced by local stringsEI. 

First of all we need to find the energy-momentum tensor of the string 
network as it evolves over time. JSpr this we adopt a model first developed by 
Albrecht, Battye and Robinson EEL The string network at any time is taken 
to be a gas of straight string segments of length L and moving with velocity 
Vrms ■ Then we evaluate the energy-momentum tensor of the network. For the 
evolution, we follow the one-scale model described earlier. A tricky part of the 
problem is to allow some of the strings to decay. This is done by eliminating 
some of the segments with time in such a way as to preserve the scaling form 
of the string energy density. Once the energy-momentum of the string network 
is obtained, we feed it into the program CMBFAST to evaluate the anisotropy 
and also the power spectrum of density fluctuations. 

Some of the results are shown in Fig. |[ (I am assuming that you have 
already been introduced to the theory behind the microwave background mea- 
surements.) Clearly there is only one Doppler peak and it occurs at I ~ 400 
and this does not agree with observations. 

It is known that larger values of Q bring down the position of the Dopnler 
peak. So Pogosian has recently investigated local strings in a closed modeltJ. 
The result that the peak position only scales as f2 -1 - 58 would seem to imply 
that we would need to go to ~ 3 in order to bring down the peak position by a 
factor of two or so. However, this result does not apply to perturbation sources 
such as cosmic strings where the anisotropy is affected by the gravitational 
potentials produced after last scattering. Simulations by Pogosian til show 
that the peak position shifts to I ~ 200 for Q ~ 1.3 (see Fig. ||). Furthermore, 
the closest fit to the data occurs for G/i = 1 — 2 x 10 -6 . While this matches 
the GUT scale, the poor fit to the data means that GUT strings cannot be 
solely responsible for producing the required density inhomogeneities. 

4-2 Gravitational wave background 

When local strings oscillate at very low frequencies, they can only radiate 
very low energy particles. Frequencies set by cosmic scales are ~ 1/t and the 



24 



Figure 2: The spectrum of microwave background fluctuations for the network of wiggly 
strings with L = 0.1. The solid and dotted lines correspond to two values of the network 
parameters. The observations with error bars are also shown. 




Sgure 3: As in Fig. ti but for a closed universe with Q = 1.3. See the article by Pogosian 
for fuller account of the parameters used and the underlying assumptions. 
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only known particles that can be radiated are massless - photons and gravitons. 
(Neutrinos are assumed to have a non-zero mass.) Generally the (local) strings 
are made of neutral fields and hence do not radiate significantly in photons fj 
Then the primary decay channel for strings is into gravitons []. 

To estimate the power emitted in gravitational radiation by a string loop 
of length L, it is most convenient to use the quadrupole formula. The only 
dimensional parameters available to us are G, [i and L. However, we know 
that the power loss into gravitational waves is proportional to one factor of G 
and two factors of /j, since the power is the square of the quadrupole moment. 
Now since the power has dimensions of mass squared, we must have 

P ~ <V (90) 

and the length of the loop drops out from the estimate. In the cosmic scenario 
we are interested in the average gravitational energy emitted per unit time by 
all strings. So we write: 

P = TGfi 2 (91) 

where T is a numerical factor. .By examining the power emitted by many 
different loops, one finds T ~ 100LL3. 

To find the present gravitational wave background amplitude and spectrum 
one needs to sum over all strings in the network, include appropriate redshift 
factors etc.. The final result isE3: 

/ 18^ 2 (/3-l) 2 ^ ( Air 

ns[Uj) - (3-/?) S in[(2-/?)7r]U/^W ( ' 

where v is a parameter that sets the amplitude for the number density of loops 
and 4/3</3<<2isa parameter that characterizes the typical spectrum of 
gravitational radiation emitted by loops, to is the present epoch. The peak of 
the spectrum is at 

47T 

Wpeafc ~ rcHTo " (93) 

The strongest constraints on the amplitude of the gravitational wave back- 
ground arise from the timing of the millisecond pulsar. In 1993 the constraint 
was: 

fl g < 4 x 10- 7 h- 2 (94) 



b If the strings are superconducting they may carry electrical charges and currents, and then 
they would lose energy by emitting photons. 

c Global string loops decay primarily by emitting Goldstone bosons. 
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The constraint gets more stringent the longer the millisecond pulsar is ob- 
served without encountering noise in its timing. Today the constraint would 
be slightly stronger than in 1993. In terms of G/x, this gives 

Gfi < 4 x 1(T 5 (95) 

The emission of gravitational radiation from strings also contributes to 
the energy density prior to nucleosynthesis. If the gravitational radiation con- 
tribution to the content of the universe exceeds about 5% , it will interfere 
unacceptably with BBN. This gives another constraint on the parameter Gfi 
that is slightly better, though comparable to, the one in eq. (p5|). 

4-3 Gravitational lensing 

Spacetime distortions due to cosmic strings would lead to gravitational lensing 
of background sources. The angular separation of images can be deduced from 
the value of the conical deficit angle to be: 

5<j) = 8ttGh-^— (96) 
a + 1 

where d is the distance between the source and the string and I is the distance 
between string and observer. This leads to an image separation of about 5 arc 
sec from GUT strings and is quite large by astronomical standards. 

Initially it was thought that a string would produce a line of lensed images 
of sources and this would be a unique signature. However, since the strings 
are wiggly, this conclusion is no longer obvious. Instead one needs to examine 
the lensing of sources due to a more realistic (wiggly) string. This was done 
by de Laix et. al. Ed and the results are shown in Fig. ||. 

4-4 Other signatures 

Cosmic strings and other topological defects have been considered as sources 
for ultra high energy cosmic rays and gamma ray bursts. 

5 Zero modes and superconducting strings 

The interaction of fermions with strings can often lead to "zero modes" i.e. 
zero energy solutions to the Dirac equation in the string background. If these 
fermions carry electric charge, electric fields along the string (equivalently, 
string motion through a magnetic field) can produce electric currents on the 
string which will persist even when the electric field is turned off. Hence such 
strings are said to be "superconducting" . 
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Figure 4: The projection of a long wiggly string on the sky is shown together with a variety 
of point sources (hatched circles). The images of the point sources are shown as unfilled 
circles. 
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The current on a superconducting string depends on the initial current on 
it at formation and also on the flux of magnetic field that it has traversed. 
The motion of the current carrying string leads to electromagnetic radiation 
that can dissipate the string energy. At the same time, the electromagnetic 
radiation can cause spectral distortions in the cosmic microwave background. 
This places constraints on the density of superconducting strings in the early 
universe and also the amount of current on them if they exist. 



6 Magnetic monopoles 

6.1 Topology: TT2 

Magnetic monopoles are formed when the vacuum manifold has incontractable 
two spheres. The relevant homotopy group is ^(G/if). As discussed in the 
string case, if ^(G) = 1 = tti(G), then 

n 2 (G/H) =n 1 (H) . (97) 

Hence we get magnetic monopoles whenever the unbroken group contains in- 
contractable paths. We know that today the unbroken group contains the U(l) 
electromagnetic gauge symmetry which contains incontractable paths. Hence, 
any GUT predicts magnetic monopoles. 



6.2 Example: SU(2) 

The 't Hooft-Polyakov monopole is based on an SU(2) model: 

L = -\x« v X a » v + \(D^ a ) 2 - V(*) , (98) 

where $ is an SU(5) adjoint scalar field, X" (a = 1, 24) are the gauge field 
strengths and the covariant derivative is defined by: 

D^ a = d^ a -ie[X^<b] a (99) 

and the group generators are normalized by Tr(T a T{,) = <5 a fc/2. Once $ acquires 
a VEV, the symmetry is broken down to U(l), consisting of SU(2) rotations 
that leave <£> invariant. 

It is possible to give a Bogomolnyi type derivation of the monopole solution 
in the case that the potential vanishes i.e. V(&) = 0. The explicit solution in 
spherical coordinates islij: 

1 Cr 

$ = ^W>T lKTa (100) 
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x j ( Cr \ 

^=e\—[l- . ur r .(<M,j = 1,2,3), (101) 
J er z \ smh(Gr) J 

This V = monopole is referred to as the Bogomolnyi- Prasad-Sommerfield 
(BPS) monopole. For V 7^ 0, the profile functions can be found numerically. 

The mass of the BPS monopole also follows from the Bogomolnyi analysis 
and is: 

M = -^m v (102) 

where my is the mass of the massive vector bosons after symmetry breaking. 
In terms of the vacuum expectation value r\ of <I>, my = er\. 



6.3 Electroweak monopoles 

The standard electroweak model does not contain the suitable topology to 
have magnetic monopole solutions. This is because iri(G) ^ 1 and the in- 
contractable loops in the unbroken U(l) are also incontractable in G. Yet, as 
originally discovered by Nambuta, the model does contain magnetic monopoles 
that are confined by strings. (This is closely analagous to the picture of quarks 
being confined by QCD strings.) 

The standard model has a doublet SU(2) gauge fields W£ and the hy- 
percharge gauge field Y^. The SU(2) and U(l) coupling constants are denoted 
by g and g' and the vacuum expectation value of $ is denoted by -q. The 
weak mixing angle W is defined via tan# u , = g' jg. The asymptotic configura- 
tion of the electroweak monopole can be written in spherical coordinates (with 
r — > oo) as: 

*=4L^U (103) 



V2 \sm(6/2)e^ 
gW£ = -e abc n b 9 Al n c + icos 2 w n a (& - 0„*t $) (104) 
g% = -isxn 2 6 w {& - $) (105) 

where 

Note that the configuration has a singularity along the negative z axis - 
the lower component of <E> is not single- valued here. This is the location of a 



real string. As described in Sec. 3.3, this string is the semilocal string when 
&w = 7r/2 and the electroweak Z-string for general 9 W . 

The electroweak monopole is not a static solution of the equations of mo- 
tion but only a field configuration. One way to see this is that to note that the 
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Table 2: The quantum numbers (rig, 713 and ni) on stable 5(7(5) monopoles are shown 
and these correspond to the SU(3), SU(2) and U(l) charges on the corresponding standard 
model fermions shown in the right-most column. 



n 


n 8 /3 


n 3 /2 


m/6 




+ 1 


1/3 


1/2 


1/6 


(u,d) L 


-2 


1/3 





-1/3 




-3 





1/2 


-1/2 




+4 


1/3 





2/3 




-6 








-1 





string is pulling on the monopole and this unbalanced force will tend to accel- 
erate the monopole and to shorten the string. Eventually the monopole will 
annihilate the antimonopole at the other end of the string. However this fact 
may not diminish the importance of the monopole - the important aspect to 
consider is the lifetime of the elcctroweak monopole and antimonopole and to 
evaluate if they can lead to some experimental signatures c3. At the moment, 
it seems unlikely that the electroweak monopole can survive long enough to be 
seen in accelerator experiments though this issue needs further exploration. 

6.4 SU(5) monopoles 

The Grand Unified symmetry breaking ST/ (5) -> [ST/ (3) x SU(2) x U(l)]/Z 6 
leads to topological magnetic_monopoles. The magnetic monopoles in this 
model have been constructed E3. Here I will not go into the details of the con- 
struction. Instead I will only point out an interesting feature of the spectrum 
of stable magnetic monopoles. 

The winding one monopole in the model has SU(3), SU(2) and U(l) 
charges. Then one can construct higher winding monopoles by assembling 
together unit winding monopoles. In some cases, the assembly will be stable 
while in other cases the higher winding monopole will be unstable to decay- 
ing into lower winding monopoles. It turns out that, for a broad range of 
parameters, the stable monopoles are the winding ±1, ±2, ±3, ±4, and ±6 
monopoles. The magnetic charges on these can be calculated. As shown in 
Table II an interesting feature is that the known quarks and lcptons have pre- 
cisely the same spectrum of charges in the electric sector! This observation has 
led to the possibility that it migb-t.be possible to understand the fundamental 
particles as magnetic monopoles 113. 
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6. 5 Cosmology of monopoles 

Once magnetic monopoles are produced at a phase transition in the early uni- 
verse, their energy density will redshift like pressureless matter p m oc 1/a 3 . 
The radiation redshifts as p 1 oc 1/a 4 . Hence magnetic monopoles will start 
dominating the universe at some epoch. The exact epoch of domination will de- 
pend on the details of monopole formation and their evolution. However, even 
if we assume one monopole per causal horizon at formation, GUT monopoles 
would have started dominating wclLhcfpre big bang nucleosynthesis, which is 
clearly in conflict with observations E3'Ej. 

Another bound on the monopole density arises from the observation that 
galaxies have magnetic fields that would accelerate any magnetic monopoles 
that are present. In this way the monopoles would dissipate the magnetic 
field. This "Parker bound" leads to a constraint that is stronger than the 
cosmological constraint for lighter monopolesEl 

The present constraints on monopoles imply that magnetic monopoles 
must be more massive than about 10 10 GeV if they are to be present in cos- 
mology. For certain types of monopoles {eg. those that catalyse proton decay), 
the constraints from stellar physics can be much stronger. 

The constraints on magnetic monopoles are clearly in conflict with the 
result that the GUT phase transition must have produced magnetic monopoles. 
Hence either the GUT philosophy is incorrect, or the standard FRW cosmology 
is incorrect. The popular solution to the monopole over-abundance problem 
is that the standard cosmology should be modified to include an inflationary 
phase of the universe. Then the monopole producing GUT phase transition 
occurs during the inflationary stage and any monopoles that were produced 
are rapidly diluted to insignificant densities. Other possible solutions are to (i) 
modify particle .physics such that the monopoles get connected by strings and 
then annihilate, E3 (ii) never have an epoch where the Grand Unified symmetry 
is restored EH, and (iii) produce domain walls that sweep away the monopoles 



7 Further reading 

The topological defects literature is vast and the citations have not done jus- 
tice to all the work that has been done. Fortunately a number of excellent 
review articles and books have been written that can be a guide to the original 
literature. Here are some of these articles: 

• "Solitons and Particles", eds. C. Rebbi and G. Soliani (World Scientific, 
Singapore, 1984). This is an invaluable resource for learning the subject 
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and for reading most of the classic papers on the field theoretic aspects. 



• "Cosmic Strings and Other Topological Defects" , A. Vilenkin and E.P.S. 
Shellard, Cambridge University Press (1994). This is a comprehensive 
description of the field. 

• M.B. Hindmarsh and T.W.B. Kibble, Rcpt. Prog. Phys. 58, 477 (1995). 
An excellent discussion of the many different facets of cosmic strings. 

• P. Goddard and D. I. Olive, Rep. Prog. Phys. 41 1357-1437 (1978). A 
classic on magnetic monopoles. 

• S. Coleman, "The magnetic monopoles, fifty years later" , in: The Unity 
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